Collective dynamics of vortices in trapped Bose-Einstein condensates 
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We have calculated collective mode spectra for three-dimensional, rotating Bose-Einstein con- 
densates in oblate harmonic traps using the microscopic Bogoliubov-deGennes field theory. For 
condensates with N v vortices, N v Kelvin-Tkachenko mode branches are obtained. The features of 
these modes are compared with those predicted by a classical point vortex model. We have created 
movies to visualize the motion of the vortices corresponding to the Kelvin-Tkachenko waves. 
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I. INTRODUCTION 

Vortical flows stand out among the most fascinating re- 
alizations of hydrodynamics. In equilibrium, an isolated 
vortex filament often settles into the shape of a straight 
line. Perturbations to such equilibrium cause helical vor- 
tex displacement waves to propagate along the length of 
the vortex filament. These excited states of vortices are 
known as Kelvin waves, named after Lord Kelvin who 
derived a dispersion relation for such vortex excitations 
PQ. Pitaevskii predicted similar vortex excitation modes 
to exist in superfluids [2]. Kelvin waves of quantized 
vortices have been observed both in superfluid helium 
[3] and in trapped atomic Bose-Einstein condensates [4] . 
However, measurement of the long-wavelength superfluid 
kelvon dispersion relation remains an experimental chal- 
lenge. 

For two or more vortex filaments, the stability of dif- 
ferent vortex configurations becomes a matter of impor- 
tance. Although rarely realized in classical hydrodynam- 
ical systems, the quantization of circulation in quantum 
systems, together with energetic considerations, leads to 
the emergence of stable vortex lattice configurations in 
quantum liquids. Considering equilibrium states of flux 
vortices in superconductors, Abrikosov predicted them 
to be arranged in a triangular lattice structure [5] , corre- 
sponding to the most efficient way of packing cylinders. 
In addition to superconductors, triangular Abrikosov vor- 
tex lattices also emerge in rotating scalar superfluids 
[7] . For vector superfluids more exotic lattice struc- 
tures appear [8HI4], 

The stability of multivortex configurations against per- 
turbations was addressed in the classical context by 
Thomson who studied the stability of 3 through 7 vor- 
tices arranged in a single orbital ring pattern [15 . Have- 
lock extended this result finding the 7 vortex case to be 
marginally stable having an excitation mode with pre- 
cisely zero frequency (modes with negative frequencies 
corresponding to instabilities), further showing that an 
arbitrary number of vortices in a ring configuration can 
be stabilized by adding sufficiently strong vortex in the 
centre of the ring [IB]. An exhaustive analysis of the 
transverse normal modes of finite vortex arrays and their 
stability was presented by Campbell [17]. Tkachcnko 
studied the stability of an infinite triangular vortex lat- 



tice in a superfluid helium and predicted the existence 
of a transverse normal mode oscillation now known as 
Tkachenko wave [15] ■ Due to the end-cap effects in ro- 
tating superfluid helium systems, the modes calculated 
by Campbell |17) did not facilitate a quantitative com- 
parison with the photographic movies of few-vortex exci- 
tation modes observed by Yarmchuck and Packard |19j . 

The experimental realization of dilute-gas Bose- 
Einstein condensates [20-22J opened up the possibility 
to controllably create and image individual quantum vor- 
tices [23] and vortex lattices p4"H27] and to study the nor- 
mal mode oscillations of such vortex arrays [7] . The vor- 
tex precession observed in the experiment by Anderson et 
al. may be viewed as fundamental Kelvin wave motion 
of a vortex in a Bose-Einstein condensate [28]. Brctin 
et al. excited a Kelvin mode of higher axial wave num- 
ber by exploiting its resonant coupling to the quadrupole 
mode of the condensate [4] . Tkachenko waves were cre- 
ated by Coddington et al. in rapidly rotating vortex lat- 
tices [55]. Smith et al. observed a gyroscopic tilting 
mode of a vortex array [50 which has one axial node and 
belongs to the common mode Kelvin-Tkachenko mode 
branch [31]. Recently, Weiler et al. [33] observed dy- 
namics of vortex dipoles and Freilich et al. applied con- 
tinuous imaging methods to obtain trajectories of such 
vortex-antivortex dipoles [33j . In three-dimensions such 
vortex dipoles may break into vortex loops via the Crow 
instability [34"r[3Tj . Similar techniques could be used to 
experimentally discover the complete set of low-lying col- 
lective Kelvin-Tkachenko mode branches of few-vortex 
arrays. 

Here we calculate the low-lying normal modes of 
oblate, three-dimensional, Bose-Einstein condensates 
hosting up to 19 vortices, using the microscopic 
Bogoliubov-de Gennes theory. This model fully accounts 
for the internal structure of the vortex cores and the cou- 
pling of vortex motion to the compressional sound waves. 
We solve the equations governing the three-dimensional 
excitation modes of such vortex arrays without forc- 
ing any symmetries in the calculations. The Kelvin- 
Tkachenko mixed modes of few-vortex single orbital ar- 
rays for three-dimensional prolate condensates were cal- 
culated in [31]. Here we extend those results to oblate 
systems with vortex lattices having up to 3 concentric 
vortex orbitals. 
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FIG. 1. (Color online) Ground states of a Bose-Einstein con- 
densate in rotating harmonic traps. Each plot is a density 
isosurface viewed along the rotation axis. The states (a)-(f), 
respectively accommodating 1, 2, 3, 7, 12, and 19 vortices, are 
being rotated at respective angular frequencies, Q' = Q,/lu± = 
0.45, 0.495, 0.497, 0.58, 0.69 and 0.79. The expectation values 
of the orbital angular momentum per particle L' — (L) / (hN) 
of these ground states are 0.97, 1.8, 2.3, 4.2, 6.2, and 9.2 
respectively. 



II. MODEL 

We employ the Bogoliubov-deGennes (BdG) wave- 
function formalism to calculate the excited states of the 
three-dimensional harmonically trapped rotating Bose- 
Einstein condensates whose multiply connected ground 
state wavefunction <fi(r) contains N v single quantum vor- 
tices. This amounts to solving the eigenvalue problem 
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where u q (r) and v q (r) are the complex valued quasi- 
particle wavefunctions labelled by q, uj q are the cor- 
responding eigenfrequencies and the chemical potential 
/j = (C(r) — g\cj)(r)\ 2 ) where £(r) is a single particle oper- 
ator defined in Eq. ([5| and g is a coupling constant. The 
numerical challenges encoded in Eq. |l]) are largely due 
to the near-singular spectrum of the BdG operator and 
the fact that for three-dimensional systems the matrix 
representation of Eq. ([!]) becomes very large. For single 
vortex lines, cylindrical symmetry is often deployed to 
reduce the computational cost of solving Eq. ([I]) [351 EH] • 
However, when more than one vortex is present in the 
system, the cylindrical symmetry of the ground state is 
lost and one has to solve the full three-dimensional prob- 
lem, which allows correct treatment of the axial structure 
of the vortex modes. We have used a finite-element dis- 
crete variable representation [40l [41] to obtain a sparse 
matrix representation for the BdG operator of this three- 
dimensional system and have used a parallelized Arnoldi 
iteration [3TJ H2] to solve the eigenvalue problem. 

The dynamics of the calculated collective modes is vi- 



sualized by tracking the evolution of the perturbed con- 
densate density \4> q (r,t)\ 2 as a function of time t, where 



^(r,*) = [0(r)+^(r ) t)]e-^/fi 



(2) 



and the perturbation caused by the excitation of the col- 
lective mode q is 



6,(r,t) =u q {r)e- l ^ t + v* q (r)e l ^ t 



(3) 



The population of the quasiparticle excitation mode is 
controlled by the number p in Eq. We have used 

p = 0.3 in this paper. 

Prior to solving Eq. (JlJ we must obtain the rotating 
ground state wavefunction 4>(r) by solving to high accu- 
racy the Gross-Pitaevskii equation 
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and the constant g = 4irh 2 a/m determines the strength 
of s— wave interactions between particles of mass m in 
terms of the scattering length a. The particles are con- 
fined by an external harmonic potential 



V cxt (r) = m (ui 2 ± r 2 + lu 2 z z 2 ) /2 



(6) 



with transverse u>± and axial w z frequencies. The sys- 
tem is rotated at an angular frequency f2, and the z- 
componcnt of the orbital angular momentum operator is 
L z = —ih(xd y — yd x ). The normalization of the wave- 
function, J v \<fi(r)\ 2 dr — N, determines the number of 
particles N in the system. Here we consider iV w 3 x 10 4 
87 Rb atoms confined in a harmonic trap with an axial 
trap frequency w 2 = 27r x 100Hz and an aspect ratio 
ui z /lj± = v°- The dimensionless interaction parameter 
gN = 1000 hujj_a,Q, where ao = {h/mwx) 1 / 2 . 



III. KELVIN-TKACHENKO MODES 

To facilitate the discussion of the collective vortex dis- 
placement modes, the corresponding ground states 4>{r) 
need to first be specified. Figure [T] shows condensate den- 
sity isosurfaces viewed along the rotation axis for six dif- 
ferent vortex systems, calculated using Eq. Q. The re- 
duced external rotation frequency O' = 0/ ui±_ and the ex- 
pectation value of the reduced orbital angular momentum 
V = (L z )/(Nh) are marked in the upper and lower right 
corners of each frame, respectively. The rotation causes 
the condensate to expand radially due to the centrifugal 
effect. Simultaneously the cloud becomes flatter in the 
axial direction. This centrifugal effect is fairly small in 
the single orbital cases (a)-(c) but becomes prominent 
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FIG. 2. (Color online) Instantaneous density isosurfaces of the 
condensate with 7 vortices, perturbed by the different lowest 
vortex displacement modes. The view is along the rotation 
axis. The excitation frequency uj' , expressed in units of the 
transverse trap frequency, of each perturbing mode is indi- 
cated in the upper right hand corner of each frame. The dis- 
crete rotation symmetries of these modes are inherited from 
the sixfold symmetry of the triangular equilibrium vortex lat- 
tice. The central vortices in (f) and (g) gyrate in the opposite 
directions [see supplementary material [SS]]. The number in 
the bottom of each frame corresponds to the mode label in 
Table I. 



for the two (d) and (e) and three orbital (f) vortex ar- 
rays. Although the superfluid circulation of the vortices 
is quantized in integer multiples of h/m, the orbital an- 
gular momentum per particle (L z )/N is not an integer 
multiple of h and hence the quantization of orbital angu- 
lar momentum is not exact despite the topological dis- 
creteness of the vortices. In contrast, for axisymmetric 
multiquantum vortices the orbital angular momentum is 
quantized in integer units |71j . 



A. General considerations 

Let us first restrict our discussion to the subset of 
modes (N v ,n z = 0), which have no axial nodes. For 
such modes all vortices remain straight and parallel to 
each other and the vortex dynamics is effectively two- 
dimensional, although the excitation frequencies do de- 
pend on the full three dimensional structure of the con- 
densate. Irrespective of the number of vortices N v > 
in the ground state, two vortex modes, which we refer to 
as (C) and (T), always appear in these systems. 

The simplest vortex displacement mode shared by all 
vortex configurations is the common mode or center-of- 
mass excitation (C) of the vortices, in which the whole 
vortex array executes retrograde circular motion, in the 
rotating frame of reference, around the trap center as a 
rigid body. Here we have deviated slightly from the ter- 
minology of Campbell [T7]. By the common mode we 
refer to the fact that the vortex motion is common to 
all vortices in a single mode, where as Campbell uses 
the term common mode to refer to fact that similar vor- 
tex mode exists for all vortex configurations irrespective 



of the number of vortices in the system. Apart from 
the obvious difference in excitation energy, the (C) mode 
may appear similar to the usual counter-rotating dipole 
mode, in which the vortex lattice stays stationary and 
the condensate density executes circular motion around 
the trap center. In contrast, in the vortex common mode 
(C), which only exists for N v ^ 0, the lattice rotates with 
respect to the trap and the stationary condensate. 

The second universal collective excitation is the trans- 
verse Tkachenko mode (T), which typically has the low- 
est excitation energy in the system and reveals itself as 
an azimuthal torsional oscillation of the vortex lattice 
[MI291 I3 HI45H57] . In this mode, the individual vortices 
travel (in the rotating frame of reference) along elliptical 
periodic orbits around their own equilibrium positions, 
their combined motion resulting in the collective mode 
of the whole vortex lattice. 

In addition to these two universal vortex displacement 
modes, for systems with at least two concentric vortex 
orbitals Campbell [17| finds "quadratic" (Q) modes in 
which the central single vortex or cluster of vortices exe- 
cutes rigid body oscillation while the vortices in the out- 
ermost orbital are able to move relative to each other. 
There are two different quadratic (Q) modes correspond- 
ing to the opposite sense of rotation of the central vortex 
cluster. Further following the notation of Campbell [T7] 
we label the modes, for which the central vortex remain 
stationary, as rational (R) modes. 

We observe that a generic collective vortex mode can 
be constructed by first grouping the vortices into sub- 
arrays and then considering the motion of the vortices in 
each sub-array as an independent collective mode. The 
vortex motion of the whole lattice is then constructed as 
a superposition of the collective modes of the sub-arrays. 
This is further illustrated in the end of the Sec. III. 

The underlying ground state vortex lattice has 6-fold 
symmetry, which is reflected by the n— fold symmetry of 
the Bogoliubov modes [55]. Furthermore, each Kelvin- 
Tkachenko branch can be uniquely labelled using a set of 
N v vortex displacement phase shifts <di for the vortices 
|31) . However this, like any means of uniquely identifying 
the individual modes, quickly becomes cumbersome for 
large numbers of vortices. 

We also note that the orbital motion of a single vortex 
around the trap centre admits three equivalent descrip- 
tions. It can equally be viewed as the long wavelength 
axially nodeless n z — Kelvin wave, the single vortex 
Tkachenko wave, or the rigid body common mode of the 
vortex lattice with only one vortex. 



B. Details of vortex modes 

Table I lists the frequencies of the calculated axially 
nodeless n z = vortex displacement modes for vor- 
tex numbers N v = 1,2,3,7,12, and 19. The Kelvin- 
Tkachenko mode branches are sorted according to in- 
creasing energy of the lowest mode in each branch and are 
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FIG. 3. (Color online) Instantaneous density isosurface of a 
condensate excited in a first axial n z = 1 Kelvin-Tkachenko 
mode corresponding to the mode branch I of the seven vor- 
tex system in Table I. The orbital vortices are symmetrically 
tilted toward the central vortex. For the full dynamics of this 
mode we refer to the attached movies [58] . 



referred to by integers. The pure Tkachenko modes (T), 
common modes (C), quadratic modes (Q) and rational 
modes (R) are labeled in the table. We have produced 
a movie featuring the collective dynamics of the vortices 
for each tabulated mode, see the supplementary material 



TABLE I. Axially nodeless n z = excitation frequencies uj q 
of vortex configurations with N v vortices. Tkachenko (T), 
common mode (C), quadratic (Q) and rational (R) modes are 
marked by the respective labels. See supplementary movies 
for the dynamics of each tabulated mode |58j . 
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The motion of N v vortices in two-dimensions can be 
approximated using a classical point vortex model |59U60j 



iT 



v ~dt 



dH 



(7) 



where z v = x v + iy v are complex functions of the spa- 
tial vortex coordinates x v and y v , Y v are the strengths 
of the vortices, and H = — J2v^ v l°g(l^ — z r/\)/^ 
is the potential driving the vortex motion when there 
are no boundaries. For confined systems more general 
form of H is used. This is a set of N v first order dif- 
ferential equations. Seeking small perturbations to the 
vortex positions leads to TV,; by matrix equation with 
N v eigenvalues corresponding to the N v normal modes of 
the vortex system. Splitting Eq. ([7]) into two real valued 
equations 



T w x v = dH/dy v 
r„y„ = ~dBjdx v , 



(8) 



yields a system with 2N V eigenmodes, only half of which 
are linearly independent the other half being related to 
the former by a simple rotation |17j . From Eqns ([8| it 
can be seen, by associating x v and y v with the canonical 
coordinate and momentum of a Hamiltonian system, that 
the phase space of this 2D vortex system corresponds to 
that of an equal number of ID massive particles. For a 
Newtonian point mass particle, the resultant force acting 
on it is proportional to the acceleration of the particle and 
therefore in two-dimensions a system of N p particles have 
2N p normal modes. However, here the Magnus force, 
which drives the motion of a vortex, is proportional to 



the velocity and therefore the vortex system has only N v 
independent low-lying normal modes. 

Equation. ^ can also be used to reproduce the vortex 
motion of the axially nodeless Kelvin-Tkachenko collec- 
tive modes by appropriately choosing the vortex poten- 
tial and the initial vortex positions. However, the result- 
ing point vortex dynamics of this chaotic Hamiltonian 
dynamical system is highly sensitive to the initial vor- 
tex locations. It therefore appears that these few-vortex 
superfluids may be particularly useful systems for stud- 
ies of quantum chaology [61] and, in particular, for the 
study of the role of wavefunction phase singularities to 
the behavior of quantum systems whose classical counter- 
parts exhibit chaos. The point vortex model was recently 
used to obtain quantitative agreement with experimen- 
tal observations for the case of a single vortex-antivortex 
pair in a harmonically trapped Bose-Einstein condensate 



1. Single orbital modes 

The Kelvin-Tkachenko modes of 1, 2, and 3 vortex sin- 
gle orbital systems are qualitatively equivalent with those 
obtained for elongated u> z /lu± = 0.2 condensates [31) . 
Briefly, the single vortex system has only one Kelvin- 
Tkachenko mode branch corresponding to the Kelvin 
waves propagating along the axis of the quantum vortex 
[64ll67j . The varicose vortex waves [68] do not involve 
the motion of the vortex phase singularities and are not 
discussed further here. The two-vortex system has two 
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FIG. 4. (Color online) Snapshots of the time evolution of the 
condensate with 19 vortices and perturbed by the rational 
mode 14 in Table I. The times r marked in the frames (a) - 
(f) are in units of 2n/uj q . 



mode branches corresponding to the in-phase (C) and 
out-of-phase (T) motions of the two vortex filaments. 
In addition to the Tkachenko and common modes, the 
three- vortex system has a mode branch with 27r/3 phase 
shift between each pairing of vortices [HO [3T] . Both the 
axially nodeless n z = and first axially excited states 
n z = 1 of the single orbital vortex systems are included 
in the supplementary videos |58| . 

Using the notation (N v , mode branch, uj q /uj±) the first 
axially excited n z = 1 Kelvin-Tkachenko modes have fre- 
quencies (1,1,1.532), (2,1,1.536), (2,2,1.658), (3,1,1.600), 
(3,2,1.690), and (3,3,1.605). The frequencies of these 
modes are approximately half of the harmonic oscilla- 
tor level spacing in the axial direction, oj z /2 = \ /r 2w±. 
By tightening the trap sufficiently in the axial direction, 
practically all axial Kelvin-Tkachenko modes can be sup- 
pressed |69j . Dynamics of the modes for the four- vortex 
system were presented in [3"T] . 



2. Multiorbital modes 



FIG. 5. (Color online) Two different decompositions of the 
same collective mode of a 19 vortex array into collective modes 
of sub-arrays, (a) six three vortex modes (b) three seven 
vortex modes. 



mode corresponding to the same branch as the mode in 
the frame (e) (cf. N v = 7, mode 1 in Table I). Here the 
central vortex remains straight while the other vortices 
are tilted performing gyroscopic motion about their re- 
spective equilibrium positions. For the characterization 
of the collective vortex modes with 12 and 19 vortices 
we refer the reader to the movies included in the sup- 
plementary material [5 8) . Figure [4] illustrates one such 
mode for the case of 19-vortex, three-orbital, configura- 
tion (N v = 19, mode 14 in Table I). Figures [5] (a) -(f) 
show snapshots of the condensate density isosurface per- 
turbed by this mode for times r = {0.15, 0.33, 0.58, 0.64, 
0.82 and 1.00}x27r/w s . The vortex motion in this mode 
can be broken down to six equivalent three-vortex clus- 
ters, see Fig.[5^a), each executing a three vortex collective 
motion (Table I, N v — 3, mode 2) while the central vor- 
tex remains stationary. Alternatively, it can be viewed 
as a superposition of three Tkachenko modes (Table I, 
N v = 7, mode 1), the 12 vortices in the outermost or- 
bital being divided into two interlaced groups and the 
inner 6 vortex orbital forming the third sub-array and 
the central vortex being shared by all three sub-arrays, 
see Fig. [5jb). 



The 7- vortex system is special. Out of all stable two or- 
bital vortex arrays it has the smallest number of vortices. 
It is also known that 7 classical vortices distributed along 
a single ring pattern is the marginally stable case such 
that any higher number of vortices will be uncondition- 
ally unstable [TH1II1]- Figure [2] (a)-(g) shows snapshots 
of the condensate perturbed by each of the N v = 7 dif- 
ferent axially nodeless n z = Kelvin-Tkachenko modes. 
The excitation frequency of these modes is marked in the 
upper right corner of each frame. Frame (a) is the com- 
mon (C) mode, (e) is the Tkachenko (T) mode, (f) and 
(g) are the two quadratic modes (Q) and (b)-(d) are the 
two rational (R) modes, see Table I. 

Figure [3] shows a density isosurface of a condensate per- 
turbed by the lowest axially excited Kelvin-Tkachenko 



IV. EXCITATION SPECTRA 

Figure [6] (a)-(f) show the frequencies of the low- lying 
modes as functions of their orbital angular momentum 
relative to the condensate ground state, for the systems 
with N v = 1, 2, 3, 7, 12, and 19 vortices, respectively. The 
Kelvin-Tkachenko modes listed in Table I are contained 
within the frequency intervals marked in each frame. The 
roton-like minima at high angular momenta seed the nu- 
cleation of additional vortices to the system when they 
become resonant with the condensate as the external ro- 
tation frequency is increased [70H72] . These "hanging" 
modes are shifted to higher angular momenta as the ra- 
dial extent of the condensate and the number of vortices 
in the system increases. The center-of-mass Kohn mode 
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FIG. 6. (Color online) Calculated low-lying excitation spectra of the respective ground states displayed in FigJT] The vertical 
intervals mark the energy range of the lowest vortex displacement modes. The center-of-mass (Kohn) mode is singled out using 
a square marker. For every frame there exists a zero mode corresponding to the condensate of atoms. In frame (d) there is an 
additional zero mode related to the broken continuous rotation symmetry of the ground state. 



always occurs at frequency u>j_ — Q and is singled out in 
the figures using square markers to differentiate from the 
Kelvin-Tkachenko waves. Each spectrum also contains a 
zero mode, which corresponds to the Bose-Einstein con- 
densate of atoms. In frame (d) there is an additional zero 
mode [73] shown also in Fig. [7]where we have plotted the 
lowest Tkachenko mode (T) marked with circles and the 
lowest common modes (C) marked with squares, as func- 
tions of the rotation frequency fl. The frequency of the 
common mode monotonically decreases with increasing 
value of tt. The frequency of the Tkachenko mode first 
steeply decreases for single orbital vortex states and then 
increases for multiorbital states. In the limit ft — > uj± the 
frequency approaches zero [221112] • The inset in Fig. [7] 
shows the rotation energy per vortex QL Z /N V as a func- 
tion of the trap rotation frequency fi. The minimum 
in this curve is correlated with the excitation frequency 
of the Tkachenko mode. Such minimum in the disper- 
sion relation for the situations where the system size is 
comparable to the intervortex spacing has not been pre- 
dicted by continuum theories, which are valid in the limit 
of large vortex lattices. The Tkachenko mode with zero 
oscillation frequency implies this system to possess a dou- 
bly degenerate rotating ground state and a finite ground 
state entropy at T = 0, unless quantum fluctuations are 
sufficient to lift this degeneracy. The stability of this zero 
mode against thermal and quantum fluctuations could be 
addressed experimentally or by using dynamical finite- 
temperature methods [73]. 

The observed two zero modes suggests the possibil- 
ity of a co-existence of two Bose-Einstein condensates in 



the system. The first one results from the spontaneous 
U(l) gauge symmetry breaking and is the Bose-Einstein 
condensate of atoms. The second zero energy quasi- 
particle mode appears due to the continuous SO (2) rota- 
tion symmetry breaking. Macroscopic population of such 
Majorana-Nambu-Goldstone boson would correspond to 
a Bose-Einstein condensate of Tkachenko waves. Fur- 
thermore, for a single vortex system it might be possible 
to create a condensate of kelvons using suitable trap- 
ping parameters. In rotating vector condensates it may 
be possible to observe further ground state degeneracies 
such as simultaneous condensates of atoms, magnons and 
kelvons. 



V. DISCUSSION 

In conclusion, we have calculated the low-lying elemen- 
tary excitation modes for rotating Bose-Einstein conden- 
sates containing up to 19 singly quantized vortex fila- 
ments, using the microscopic Bogoliubov-deGennes the- 
ory. For N v vortices in the system, the spectrum con- 
tains N v low-lying Kelvin-Tkachenko collective mode 
branches. We have created movies of these excitation 
modes, showing the dynamics of the vortices when these 
collective modes arc populated. In agreement with the 
classical point vortex model, for every vortex configura- 
tion we find two universal modes. These are the fun- 
damental Tkachenko mode and a vortex lattice center- 
of-mass common mode. For configurations with more 
than one concentric orbital of vortices, we also find two 
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FIG. 7. (Color online) Excitation frequencies ui q of the low- 
est Tkachenko modes (filled circles) and the lowest common 
modes (squares) as functions of the external trap rotation fre- 
quency Q,. The inset shows the energy per vortex associated 
with the rotation of the trap as a function of the trap rotation 
frequency. 



quadratic modes for which the vortices inside the outer- 
most orbital move together as a rigid body. In addition, 
several rational modes, for which the central vortex re- 
mains stationary, are found. 

These Kelvin-Tkachenko modes could be experimen- 
tally observed using the existing experimental cold atom 
technologies. An interesting future direction is to in- 
vestigate to what extent the physics of two-dimensional 
quantum turbulence can be understood in terms of non- 
adiabatic collective motion of vortices [75] . For the em- 
ployed parameters the fundamental Tkachenko mode of 
the system with seven vortices has a zero oscillation fre- 
quency |73j . Such zero mode may be interpreted to signal 
Bose-Einstein condensation of Tkachenko wave quanta in 
the system. For other studied vortex configurations this 
Majorana-Nambu-Goldstone mode is found to have finite 
excitation frequency, lifting the degeneracy observed for 
the seven-vortex ground state. Further studies of the 
spectra of slowly rotating Bose-Einstein condensates are 
needed to investigate the robustness of the observed min- 
imum in the excitation frequency of the lowest Kelvin- 
Tkachenko mode. 
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